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We present a generalized Gross-Pitaevskii equation that describes also the dissipative dynamics 
of a trapped partially Bose condensed gas. It takes the form of a complex nonlinear Schrodinger 
equation with noise. We consider an approximation to this Langevin field equation that preserves 
the correct equilibrium for both the condensed and the noncondensed parts of the gas. We then 
use this formalism to describe the reversible formation of a one-dimensional Bose condensate, and 
compare with recent experiments. In addition, we determine the frequencies and the damping of 
collective modes in this case. 



The observation of Bose-Einstein condensation in ul- 
tracold trapped atomic vapors fl]-|| has offered the possi- 
bility to study the equilibrium and nonequilibrium prop- 
erties of these degenerate gases experimentally and to 
compare the results with ab initio calculations. The lat- 
ter have as an input only the mass and scattering length 
of the particular atom of interest, and the parameters 
involved in the experimental setup. To describe the var- 
ious equilibrium and nonequilibrium properties like for 
example topological excitations, relaxation rates, mode 
frequencies, damping rates and density profiles, theories 
have been developed at different levels of sophistication. 

At the most elementary level, the Gross-Pitaevskii 
equation already captures many of the experimentally 
observed phenomena Q . It describes in Hartree approx- 
imation the zero-temperature dynamics of the conden- 
sate, and has been used to explain and predict many 
features of these Bose-condensed systems. At the next 
level, the static and dynamic properties of the noncon- 
densed or thermal part of the gas are to be included. This 
can in first instance be done by including into the Gross- 
Pitaevskii equation the Hartree- Fock interaction with the 
thermal cloud and coupling it to an equation for the dy- 
namics of the thermal part of the gas. The latter is in 
good approximation given by a Boltzmann equation for 
the single-particle distribution function . 

All these theories describe in essence only the average 
dynamics of the gas. Near the critical region however, 
fluctuations in the order parameter are generally much 
larger than the average value of the order parameter it- 
self. Therefore, in this region, it is necessary to include 
fluctuations into a description of the trapped gas. Also 
far below the critical temperature, it can be essential in 
some cases to include fluctuations into a description of 
the dynamics of the gas. For example, to understand 
the phenomenon of phase diffusion J(| , one needs to con- 
sider fluctuations that disturb the phase of the conden- 
sate. From a fundamental point of view, the desired dis- 
sipative generalization of the Gross-Pitaevskii equation 
should obey the so-called fluctuation-dissipation theorem 



0. This guarantees that both the condensed as well as 
the noncondensed components of the gas will relax to 
thermal equilibrium. As a result, the order parameter 
fluctuates around its mean value, and the central quan- 
tity describing the dynamics of the condensate is not this 
mean value, but the actual probability distribution of the 
order parameter. With all this in mind, a unified theory 
describing the coherent dynamics of the condensate wave 
function, the incoherent scattering between the various 
components of the g well as the fluctuations around 
the average value of the order parameter has been devel- 
oped by one of us $M . 

The purpose of this Letter is two-fold. First, we show 
how the coupled dynamics of the thermal cloud and the 
condensate can be solved in a selfconsistent way, that on 
the one hand leads to the correct equilibrium distribu- 
tion of the trapped gas, and at the same time takes into 
account both mean-field effects as well as fluctuations. 
Second, we show that fluctuations can be of crucial im- 
portance when trying to understand recent experimental 
results. 

In general, a theoretical description of a trapped in- 
teracting Bose gas is possible in terms of the Langevin 
equation 
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Here Ti is Planck's constant, m is the mass of the atom, 
T4xt(x) is the external trapping potential, and T 2B = 
iwh a/m is the s-wave approximation to the two-body 
scattering matrix, with a the scattering length. This 
Langevin equation can be derived using a field-theoretic 
formulation of the Kcldysh formalism fl, and describes 
the fluctuations as well as the mean-field effects of both 
the condensed and the noncondensed parts of the gas. 
The derivation requires that the high-energy part of the 
system is sufficiently close to equilibrium that it can be 
described as having a temperature T and a chemical po- 
tential jii, and can play the role of a 'heat bath'. This 
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requirement is usually well satisfied and enters through 
the explicit expression for the imaginary part in our gen- 
eralized Gross-Pitaevskii equation, which, if the gas is 
sufficiently close to or below the critical temperature, can 
be approximated by 
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Here (3 is l/fc^T, with kg Boltzmann's constant. 

If i?(x, t) has this particular form, the trapped gas will 
relax to equilibrium, because it enforces the fluctuation- 
dissipation theorem. Note that Eq. (||) is always valid for 
the energy levels below the chemical potential, i.e., for the 
condensate, but causes the energy distribution function 
for the noncondensed cloud to relax to its 'classical' value 
N(e) — \fi(e — Therefore, it cannot describe the 

exponential decay of the density of noncondensed atoms 
at the edges of the thermal cloud. Instead, the density of 
thermal atoms decays algebraically. In addition, the clas- 
sical approximation overestimates the average number of 
atoms for the eigenstates above the chemical potential. 



Both defects can be cured |10|, but are unimportant for 
the condensate and the low-energy part of the thermal 
cloud, where most of the atoms reside, and for which our 
theory is intended to be valid. 

Finally, the correlations of the Gaussian noise 7?(x, t) 
in Eq. (|lj) are given by 

(r 7 *(x,f)r 7 (x',t')) = \v K (x)8{t - t')S(x - x') , (3) 

where the average denotes an average over the differ- 
ent realizations of the noise. The general expressions for 
i?(x, t) and T, K (x) can be found in Ref. 0, and lead to 
a Keldysh selfenergy T, K (x) equal to 
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with Ni = N(ei) the Bose distribution for the eliminated 
part of the gas, and a = p 2 /2m + V r cx t(x) — /i. It de- 
termines the strength of the fluctuations through Eq. (||) 
and is related to the damping by means of the fluctuation- 
dissipation theorem. 

At this point, we explain briefly the experimental setup 
we are considering in the rest of this Letter. It is in- 
spired by a recent experiment by Stamper-Kurn et. al 
jnj , and the ideal realization of the conditions mentioned 
above. In the experiment we are considering, a Bose gas 
is trapped and cooled to a temperature above the tran- 
sition temperature. Subsequently, a dimple is created in 
the external trapping potential, say along the z-axis, by 
means of optical techniques. This dimple is steep enough 



such that there is only one energy level in the poten- 
tial perpendicular to the z-axis. We assume the dimple 
to be well approximated by a harmonic potential with 
trapping frequency uj±. Factorizing the wave function 
<&o(x) = &(x, y)$(z), the atoms trapped inside this dim- 
ple form effectively a one-dimensional gas, with an in- 
teraction strength g — T 2B /2tt1 2 l . Here, l± — y/h/mcu± 
is the harmonic oscillator length, and <f"o the harmonic 
groundstate in the dimple. By changing the depth of 
the dimple, its lowest energy level can become lower 
than the chemical potential of the noncondensed three- 
dimensional gas. If this situation occurs, the atoms will 
condense into this ground state. Notice that during this 
process, the noncondensed gas in the three dimensional 
trapping potential will remain close to equilibrium, and 
represents the 'heat bath'. 

We now turn to a numerical solution of Eq. (|l|) 
under these conditions, using a combination of well- 
known techniques |D1 |. To ensure particle number 
conservation in the absence of the term iR(x,t), we 
use an implicit method that represents the time evo- 
lution operator exp(—iH6t/Ti) as (I + iHSt/2?i)~ 1 x 



(1 - iHSt/2h) 
fi — iR(x, t) + 



Here, 



H = - 7rV72 m + Fext(x) - 

g|$(x, t)| 2 , with <I>(x, ti) the selfconsis- 
tent average [$(x, t{) + $(x, ti + St)]/2. The numeri- 
cal method for solving Eq. (Q) can now be found from 
its solution $(x, U + St) = exp(-iH5t/h)[^(x,ti) - 
(i/h)exp(-iHU/h) ft* +St dt' exp(iHt' /H)r)(x,t')} by in- 
troducing a new noisy variable £;( x ) = exip(—iHti/h) 
L' * dt' exp(iHt' /h)r](x,t')]. The correlations of £i(x) 

are(e*(x)0(x')) = ( x )MW(x-x') + 0(<5i 2 ). The 

spatial discretization is straightforward. 

As a first application, we show in Fig. I that in the case 
of a noninteracting vapor, the density of the harmonically 
trapped one-dimensional gas relaxes to the correct equi- 
librium given by n(z) = J2 a \4>a(z)\ 2 / f3(e a — /i), where 
e a are the energy eigenvalues and <p a (z) the correspond- 
ing eigenfunctions of the Hamiltonian. The equilibrium 
is shown for several sizes of the spatial mesh, and we 
see that the density distribution converges towards the 
continuum limit given by [2l 2 z f3Tiuj z ]~ 1 / 2 [f3(y e y L t — ^)]~ 1 / 2 
where l z = ^Jh/mu) z is the harmonic oscillator length 
along the z-direction, only for rather small mesh sizes. 
This is caused by the relatively large contribution of high- 
energy states due to the classical behavior of the thermal 
cloud. We emphasize, however, that given a certain mesh 
size, the correct equilibrium corresponding to that partic- 
ular discretization is reproduced numerically. Moreover, 
it is important to realize that if we do not include the 
noise, the density of the gas would be zero. Hence, it 
is clear that the fluctuation-dissipation theorem ensures 
that the noise and the imaginary term in Eq. (Q) cooper- 
ate in order to occupy the energy levels thermally. Thus, 
including fluctuations is crucial in treating the effects of 
the thermal cloud. 
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FIG. 1. Stationary density profile of a noninteracting 
gas above the critical temperature. The time step used is 
At = 0.002 wj 1 , the discretization Ax used for the differ- 
ent lines is [1]: Ax = 0.05 l z , [2]: Ax = 0.2 l z , and [3]: 
Ax = 0.8 l z . Also shown are the predicted classical density 
profiles for curves [1] , [2] , and [3] (dashed lines) , and the con- 
tinuum result (solid line). The chemical potential of the three 
dimensionally trapped gas is /i = —30 hu) z , and its tempera- 
ture T = 400 nK. The trapping frequencies are oj z — 2tvx 13Hz 
and lo y — 2-7T x 20Hz JnJ|. The trapping frequency of the dim- 
ple is uj±_ — 2n x 500Hz, and the ground state eigenvalue 
eo = — 25 Tiuj z . The atom considered is sodium, with a scat- 
tering length a « 2.75 nm. 

Second, we calculate the damping of the breathing 
mode in a trapped one-dimensional Bose gas below the 
critical temperature. Shown in Fig. 2 is g times the value 
of the density n(0) in the center of the one-dimensional 
trap. First, the gas relaxes towards equilibrium at an ef- 
fective chemical potential \i — eo = 30 fud z , where eo is the 
eigenvalue of &o(x, y) in the dimple. The fact that the 
value of gn(Q) does not relax to 30 tko z can be explained 
by realizing that mean-field effects are included in our 
formalism, because $>(z,t) describes both the condensed 
and the noncondensed parts of the one-dimensional gas. 
This implies that the chemical potential for the con- 
densed part of the gas is effectively lowered by the in- 
teraction with the noncondensed part of the gas. After 
allowing the gas to equilibrate, at t = 0.61 s the trapping 
frequency is instantaneously changed to u>' z = 0.95 u> z , 
and the gas starts to oscillate. The frequency ui and 
the damping rate 7 of the oscillation are found by fit- 
ting to n(0,t) = n cq (0) + Sn(0) exp(— 7*)[1 — cos(wt)], 
which describes both the relaxation of the density to its 
new equilibrium, as well as the excitation and damping 
of the oscillation. The results are ui = 136.84 s _1 , which 
is close to the result of ^/5/2cj z = 129.1 s _1 expected 
from the Gross-Pitaevskii equation (L4|, 7 = 12.3 s -1 , 
g6n(0) = 0.84 huj z , and gn eq (0) = 28 Huj z . Note that the 
damping is in principle caused both by the collisions with 
the reservoir of thermal atoms in the three-dimensional 



trapping potential, which is described by E^(x), as well 
as by the nonlinear term in the Langevin equation which 
induces damping upon averaging over the different real- 
izations of the noise. The latter includes collisional, as 
well as Landau damping. 

In Fig. 3, a snapshot of the density profile is shown 
at several times during the initial growth of the conden- 
sate. These are taken from the simulation presented in 
Fig. 2. As expected, the condensate shows up as a peak 
in the density profile. Also shown are the Thomas-Fermi 
solution \p — V cx t{z)]/g for the condensate density and 
the noninteracting result for the noncondensed part of 
the gas which diverges at the critical point where the ef- 
fective chemical potential becomes equal to zero. This 
is due to the fact that in one dimension, mean-field the- 
ory completely fails near the critical temperature. Fig. 3 
shows that far enough from the condensate, the noninter- 
acting result is obtained. Again, in the center of the trap 
one obtains the Thomas-Fermi result only approximately 
due to interactions. 
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FIG. 2. Oscillations of the trapped partially 

Bose-condensed gas, and their damping, after a small and in- 
stantaneous change of the trapping frequency tj z at t = 0.61 s. 
Here, we used At = 0.002 lu' 1 , Ax = 0.05 l x , e Q = ~QQhuj z , 
and u>± — 2n x 1200Hz. The remaining parameters are the 
same as in Fig. 1. 

Finally, we consider the reversible formation of a one- 
dimensional Bose condensate, similar in spirit to the ex- 
periments by Stamper-Kurn et. al jnj. The dimple per- 
turbing the external trapping potential is now oscillating 
in such a way, that the lowest energy level in the dimple 
crosses the chemical potential of the three dimensionally 
trapped gas several times according to eo = + sin(cjt), 
with u) — 2ir s _1 and /i = — 30huj z . The calculations have 
been done with and without noise. The results show that 
without noise, the condensate evaporation and growth 
cannot be described properly. Indeed, in that case, our 
findings depend strongly on the initial conditions. More- 
over, the periodic growth occurs with decreasing ampli- 
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tude. 




FIG. 3. This figure shows the growth of the density profile 
of the trapped partially Bose-condensed gas for the simulation 
of Fig. 2. The snapshots are taken at t — 2 lo~ x , t = 4 lo' 1 , 
t = 10 uj- 1 , t = 20 to' 1 , t = 30 wj 1 , and t = 40 ojJ 1 . The 
dashed lines show the equilibrium mean-field results in and 
outside the condensate. 

To describe the growth cycles correctly, one needs to 
include fluctuations into the generalized Gross-Pitaevskii 
equation, which is to be expected since we are at times in 
the critical region. Notice that we quantitatively repro- 
duce the lagging behind of the condensate as observed 
in an experiment performed with an essentially three- 
dimensional dimple. We observed a lagging behind of 
roughly 0.1s, whereas in experiment 0.07s was measured 
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FIG. 4. Reversible formation of a condensate (solid line) 
because of a sinusoidal time dependence of the chemical po- 
tential in the dimple (dashed line). The time step used 
was At = 0.003 to' 1 , the discretization Aa; = 0.05 l z , and 
uj±_ = 2n x 1200Hz. The other parameters used are the same 
as in Fig. 1. The two dot-dashed lines correspond to simula- 
tions without noise, with 10 (top) and 0.1 (bottom) particles 
in the initial harmonic oscillator ground state. 



In conclusion, we have shown that our generalized 
stochastic Gross-Pitaevskii equation consistently takes 
into account the dissipative dynamics of a trapped par- 
tially Bose-condensed gas, and describes the equilibrium 
density profile, condensate growth, coherent dynamics, 
and damping in a unified way. In our opinion, the one- 
dimensional experiment considered here would be ideal 
for a detailed comparison between theory and experiment 
in the problem of condensate growth, because the three- 
dimensional cloud remains in equilibrium. In addition, a 
numerical study of the two-dimensional case would be of 
interest and might be used to investigate for example the 
dissipative dynamics and the formation of vortices in a 
rotating Bose gas. 
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